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Abstract 

In this work, we adopt the simplest model that spontaneously breaks super symmetry, namely, the 
minimal O'Raifeartaigh model. The effective potential is computed in the framework of the linear 
delta expansion (LDE) approach up to the order 5 2 , conjugated with superspace and supergraph 
techniques. The latter can be duly mastered even if supersymmetry is no longer exact and the 
efficacy of the superfield approach in connection with the LDE procedure is confirmed according 
to our investigation. That opens up a way for a semi-nonperturbative superspace computation 
which allows us to deal with spontaneously broken super symmetric models and encourages us to 
go further and apply this treatment to the Minimal Supersymmetric Standard Model (MSSM) 
precision tests. 
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I. INTRODUCTION 



The thrilling times of the LHC Physics we are living in opens up a great deal of issues 
connected to fundamental mechanisms and theories, in special, supersymmetry (SUSY) and 
its possible breaking mechanisms [l-lQ]. Considering fundamental principles of quantum 
fields theory, SUSY seems to be fairly well motivated as a very fundamental symmetry of 
the high-energy regime. At our accessible energies, it does not show up; it has to be broken 
at some scale much above our reachable energies, and its possible evidences at accelerator 
energies must be communicated by means of some mechanism connecting this (higher energy) 
breaking scale to our low-energy world. 

The present and the near-future LHC outcomes are crucial for the interplay between 
SUSY and the Standard Model parameters. The focus is not on SUSY itself, once we 
understand that SUSY is very likely to show up at very high-energy scales; the actual matter 
with SUSY relies on its possible breaking mechanisms and the ways the latter are driven 
down to the cutoff region of the Standard Model, namely, the TeV scale. In this framework, 
the quest for possible new SUSY violation mechanisms and a broader exploitation of the 
already known models to breakdown the fermion/boson symmetry are self justifiable 11]. 

On the other hand, LHC Physics is also refining the precision tests and the level of 
accuracy of the measurements of the Standard Model's parameters. Since perturbative 
quantum field-theoretic calculations are the way we get the phenomenological results for 
Standard Model processes to be compared with experiments, we are face to face with the 
need to go further in perturbation theory, so as to incorporate higher order corrections into 
the calculation of physical processes. 

Placed in this scenario, we are motivated to reassess SUSY breaking models by computing 
higher order corrections to their corresponding effective potentials, so as to probe the effects 
of SUSY breaking in connection with the improvement of precision tests at the LHC. It is 
clear that LHC is a collider for new discoveries rather than a precision machine; but, anyhow, 
it increases the level of the precision tests of the previous LEP. It becomes a mandatory task 
to ascertain how much loop corrections affect the pattern of SUSY breaking once we start 
off from a violation that takes place at the classical level. 

In connection with the discussion of SUSY breaking to account for the splitting of the 
masses of supersymmetric partners, we would like to point out that, very recently, a new 
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proposal of a model based on SUSY has been proposed by Alvarez, Valenzuela and Zanelli 



121 ], in which a structure of partners do not show up, although SUSY is locally realised. So, 
there is no need of SUSY breaking to split masses and the (fermionic) matter fields acquire 
mass through their coupling with some background geometry. 



In a series of previous works 13j, [l4|, we have adopted the minimal O'Raifeartaigh model 
[l| , which realises SUSY violation by means of the so-called F-terms, and we have deviced a 
technique to approach the problem with the use of superfield and supergraph techniques. To 
get a richer perturbative series, we have chosen the so-called linear delta expansion (LDE) 1 
and we have coupled this method to our superfield methods. The outcome was encouraging 
and, in view of the efficacy of the conjugating supergraph techniques with the LDE, here 
we propose to carry out a computation of the effective potential up to 0(5 2 ). Owing to 
a particularity of the LDE, calculations at this order require to take into account vacuum 
diagrams up to two loops. 

The LDE is a nonperturbative method which automatically resums large classes of terms 
in a self-consistent way, to avoid possible dangerous overcounting of diagrams. This is 
achieved by combining perturbation theory with an optimization procedure. It has a long 
history of successful applications, describing phenomenological models using quantum field 
theory at zero temperature and under extreme conditions. It has been shown that the LDE 



results go beyond the standard mean-field or large- iV approximation 
finite-iV effects. Some very interesting results can be founded in 



jy explicitly including 
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23) and in references 



therein, and strong signals of the convergence of the method can be found in [24], 125 ]. 

Usually, when two (or more) loops are present in the perturbation series, we need to 
implement a numerical calculation to perform the optimization procedure. In this case, it 
was shown in all applications cited above that the numerical results of the LDE can go 
beyond the usual resumation methods. Here we further develop the superspace applications 
of the LDE, by taking into account 0(5 2 ) terms in the effective potential expansion and 
solving numerically the optimization procedure. In particular we study the convergence of 
the method, where we contrast the numerical results with our analytical results obtained at 
0{5 1 ) using two different optimization procedures. 



1 in recent studies in the literature some authors have called the LDE method as optimized perturbation 
theory (OPT), since the method is not only an expansion in the parameter 5, there is an optimization 
procedure on the method. 
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The general structure of our paper is as follows: in Section II, we briefly review the 
application of the LDE to supersymmetric field theories, while working in superspace. In 
Section III, we employ supergraph techniques to compute the one- and two-loop diagrams 
that contribute to the 0(5 2 ) to the effective potential. All the perturbative calculations 
of Section III and the numerical results we work out are collected in Section IV. Finally, 
our Concluding Remarks are cast in Section V. An Appendix follows, where we list all the 
superspace integrals of the supergraphs evaluated in Section III. 



II. CATCHING-UP OF SUPERSPACE LINEAR DELTA EXPANSION 



The purpose of this section is a warming-up with a general presentation of the linear 
delta expansion (LDE) in the frame of (matter) supersymmetric field theories. We adopt a 
superfield approach and follow references jl^, 14, 26 1. Building up our superspace action in 



terms of chiral and antichiral supermultiplets, we start off from what we call the interpolated 
Lagrangian C 5 : 

C 5 = 8£(fi,p) + (l-8)£ (fi,p) , (1) 

where 8 is an arbitrary parameter, Co{p,p) is the free sector of the Lagrangian, and /i, p 
are mass parameters. Notice that, whenever 8 = 1, we recover the original Lagrangian. 
The 8 parameter appears in connection with the interaction terms and is so chosen to be 
the perturbative expansion parameter; this means that we do not perturbatively expand in 
terms of the coupling constant itself. The mass parameters appear in £ an d 8Cq. The 
(/!, p) dependence of Co is summed up into the propagators, whereas 8Cq is regarded as an 
insertion and is taken as a quadratic interaction. 

Let us now state our methodology. We carry out a usual perturbative expansion in 8 and, 
at the very end of the calculations, we take 8—1. At this stage of our approach, ordinary 
perturbation theory is applied and a finite number of Feynman diagrams is calculated; the 
results are essentially perturbative. However, quantities evaluated at a finite order in 8 
explicitly depend on the parameters \i and p. Therefore, it is necessary to fix them up. To 
do this, we adopt the principle of minimal sensitivity (PMS) j^]. In this framework, the 
effective potential V^j(fi, p), perturbatively evaluated to order 8 k , must be taken at a point 
where it is less sensitive to the parameters p, p. Invoking the PMS, p — Ha and p = po 
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appear as solutions to the equations 




M=M0,<5=1 



M=M>,5=1 



= 0, 



(2) 



and will come out as functions of the original coupling and fields. We then insert fiQ, Jxq in 
the expression for our effective potential V^, obtaining a nonperturbative result, once our 
propagators depend on fi, p,. 

Our method is based upon the calculation of all the diagrams up to a given order in 
5, including the vacuum diagrams. In ordinary quantum field theory, for the calculation 
of the effective potential, we do not in general worry about vacuum diagrams, since they 
do not depend on the fields. In our approach, the vacuum super-diagrams do depend on 
/j, and fi and become important to the LDE, once the arbitrary mass parameters will now 
depend on fields by virtue of our optimization procedure. Thus, in the LDE, an order-by- 
order calculation of the vacuum diagrams becomes mandatory. On the other hand, vacuum 
diagrams in superspace are identically zero, due to the Berezin integrals. In view of this 
undesirable feature, we take, from the onset, the parameters /x, p, as constant superfields and 
we keep the vacuum superdiagrams till the end of the optimization procedure. To render the 
procedure clearer, let us give below the expression for the superfield generating functional 
in the presence of the sources J and J (chiral and antichiral, respectively): 



of the propagator and, in addition to the original interaction terms, one has new bilinear 
chiral and antichiral interaction terms proportional to Sfj, and 5p. We are then lead to the 
following expression for the superfield effective action: 




with m being the original mass, M = m + ft and M — m + ji. G^ M ' M ^ is the matrix form 




where sDet(G^ M ' M ^) is the superdeterminant of G^ M ' M \ It is in general equal to one, but here 
we keep it, since G^ M ' M ^ depends on fi and p. Besides, in view of the fx and p dependence, the 
generating functional of the vacuum diagrams, Z[0, 0], is not identically one. We can define 
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a normalized functional generator as Z 



N 



Z[J,J] 
Z[0,0] ' 



and set the effective action as written in 



the expression below: 



$] = -- \n[sDet(G)} - i In Z[J , J ] + r N [$, $] 



where the sources Jq and Jq are defined by the equations 



5W[J, J] 



5J(z) lJ ~ Jo 5J{z) U=J0 6J(z 

In the first two terms (usually equal to zero) stand for the vacuum diagrams and Fjv[$, $] 
is the ordinary contribution to the effective action. 

Let us now work out the interpolated Lagrangian and the new super-Feynman rules for 
the O'Raifeartaigh model. The minimal O'Raifeartaigh model is specified by the Lagrangian: 



8W[J, J] 
5J{z) |j=Jo 



5Z[J, J] 



\J=Jo 



5Z[J, J] 
5J{z) |j=Jo 







(5) 



(6) 



C 



J d 2 9 (£$ + m$ 1 $ 2 + #$ $i) + /i.e. 



(7) 



where i — 0, 1, 2. 

Following the work of Ref. [13], to account for the nonperturbative contributions of all 
fields in the model, we apply the LDE with the matrix mass parameters and fiij. By 
adding and subtracting these mass terms to a genral O'Raifeartaigh-type action, we obtain 



£(//, fl) = £ (/i, /j) + CtatQi, (i) 



with 



d 2 9 + -MiJb&j + h.c 



(8) 

(9) 
(10) 



Coin, ft) = J 

Cint{^V) = - J d 2 6 (J^g ijk $ i $ j $ k -^iii j $ i $ J J +h.c. 

where My = m^- + and i, j, k = 0, 1, 2 are symmetrized indices. 

We now cast the superfield expansions for the arbitrary mass parameters as follows: 

My = ^ijk^k = \jk{Pk + # 2 Xfc) = KjkPk + \jkXk0 1 = Pi? + ^'^ 2 , (H) 



so that 



Mi. 



"\? + fJ>ij = {rriij + Pij) + kjO 2 = a»ij + kft 2 



(12) 



The interpolated Lagrangian ([I]) takes the form 



int J 



(13) 



where the free and the interaction Lagrangians are given by: 



r 6 



d 2 9 ( ;,«!•, + -tty^j + ) h.c. 



int 



5 5 
d2 ° ( y9ijk®i®j®k ~ ^ij^i^j 



h.c. 



(14) 



(15) 



Notice that the interaction Lagrangian now displays soft SUSY-breaking terms propor- 
tional to the p components. We treat these terms perturbatively in 8 as for normal interac- 
tion terms. 

Now, to recover the so-called minimal O'Raifeartaigh model, built up in terms of three 
superfields, with 5 = 1 (ED), we choose: 



£o = £; 

M i = a 01 = poi = a ; 
Mn = b n 6 2 = be 2 ; 

M 12 = a\2 = "2i2 + Pu = m + p = M ; 
I, #011 = 9 , 

and all other £j and Mij set to zero. With these choices, we obtain 



(16) 



C 



hit 



(17) 



This O'Raifeartaigh model possesses an (Abelian) R symmetry. The R charges assign- 
ments of the chiral superfields $o> ^i, &2 are respectively Rq = 2, R\ = and i?2 = 2. To 
keep the R symmetry of the interpolated Lagrangian, the R charges of the parameters a and 
b are R a = and Rb = 0; they are to be left the same after the optimization procedure. 

The new set of modified propagators can be read off from the free Lagrangian, which has 
also explicit dependence on 9 and 8 through p and p. Using the techniques developed in 
, the new superfield propagators read as follows: 



($ *o) 
($o$i) 
($0*2) 



(k 2 + \M\ 2 )A{k)5{ 2 + \a\ 2 \b\ 2 B(k)9 2 9 2 8 

abC(k)±D 2 D 2 e 2 5t 2 ; 
lo 

-MaA(k)5\ 2 + Ma\b\ 2 B(k)el9 2 5 



12 > 



12 ) 



E(A;)4 + |6| 2 S(fc)-D^ J D?4; 



($1$ 2 ) 


= -MbF(k)6{5{ 2 ; 


($ 2 $ 2 ) 


= (k + \a\ 2 )A(k)5f 2 + \M\ \b\ B(k)0(6(5f 2 




- -\n\ 2 br(k) 1 D 2 f) 2 5 4 ■ 




- «4ft^ n 2 A 4 — «i/ii 2 Ri't^/9 2 ^ 2 n 2 A 4 ■ 
— a^Acj-i^o^ a\o\ tf{K)-v 1 v 1 u 1 o 12 , 


\^0^2/ 


- — MnbC(k'\-D 2 f) 2 5 4: ■ 

— 1V1 VbVKj yrv J ^ 1 1 12 ' 




= ^ W-U^d^ , 


($1$ 2 ) 


= MA{k)\Dl5t 2 -m 2 B{k)\Dl9l9l5{ 2 


($ 2 $ 2 ) 


= -\M\ 2 bC{k)\Dlel5t 2 , 



where 



(18) 



A{k) ~ P(P + |M| 2 + |a| 2 ) ' (19) 



B(k) 
C{k) 
E{k) 



1 



(k 2 + \M\ 2 + \a\ 2 ) [{k 2 + |M| 2 + \a\ 2 Y - \b\ 2 ] 



k 2 [(k 2 + \M\ 2 + \a\ 2 ) 2 - \b\ 2 ] ' 
1 



k 2 + \M\ 2 + \a\ 2 



F(k) = = . (20) 

V ; (P + |M| 2 + |a| 2 ) 2 -|6| 2 V ; 

It is noteworthy to highlight the nontrivial dependence of the propagators on the param- 
eters a and b: the optimized parameters are present at the poles of the propagators, as if 
they were typical mass terms. 

The new super-Feynman rules for the vertices are: 





vertex : 


n 9 j 


d A 9; 




vertex : 




[ d 4 9; 




vertex : 


—5a 


[ dH- 


<3>i$i 


vertex : 


5b 


1 d 4 ee 2 



(21) 

We are now ready to start calculating the perturbative effective potential in powers of 5. 
To do that, we use the vertex functions defined in the expansion of the effective action and 
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'1 4^1 Sf2 

FIG. 1: Effective potential up to the order 8 1 



consider the compatible vacuum diagrams. In Ref. 13J, after the optimization procedure, it 
is seen that the order-5° contribution accounts for the sum of all one-loop diagrams. It was 
then possible to derive analytical solutions to the optimization procedure before calculating 
the superspace and loop momentum integrals. Nevertheless, in second order, we are not able 
to exhibit analytical solutions; so, we must renormalize the theory before the application of 
the optimization procedure. To do that, a numeric calculation is required. 

In Fig. [TJ we depict the diagrammatic sum of the effective potential up to the order 5 1 

Notice that, due to the ^-dependent propagators, tadpole diagrams do not vanish any 
longer, as it usually happens in superspace. The first diagram is of order 5° and we can see 
that it corresponds to the first term of the effective action, as given in (jSJ). The third graph 
is a vacuum-type superdiagram with a quadratic insertion that stems from the 9 2 component 
of the a expansion. 

n 

With our super- Feynman rules and the results of [13J, the effective potential up to the 
order S 1 reads as given below: 



i=i 

= \j d% 2 5t 2 Tr\n[P T K]5t 2 + 6 j 0^F(k) {-2gb j d 2 6% + \ \b\ 2 + h.c.} 
+S J (0 S(fc) { 4 ^ 16 ' 2 / d2de2 ^-\ a \ 2 \ b \ 2 -pM\b\ 2 + h.c.} . (22) 



In Ref. 14|, the renormalization of the effective potential up to the order 5 2 is discussed. 
At the order 5 1 , the counterterm below is needed: 

25gb f j4/1 ^ 2 ^ , ; 25gb 

K€ 

where $ _r denotes the renormalized superfield. The chiral potential is therefore renor- 



J d 4 66 2 $ R + h.c. = ^ j d 2 9% R + h.c. , (23) 



malized. The counterterm that renormalizes it depends on b, which is a solution to the 



optimization procedure. In 



141 ]. it was shown that, after the optimization, only the Kahler 



potential is actually renormalized. 
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The renormalized effective potential up to the order S 1 is 



V. 



(i) 

eff 



(4tt) 



(M 2 + a 2 ) 2 In 



1 - 



+-(M 2 + a 2 )ln 



(M 2 + a 2 Y\ 
M 2 + a 2 + 6 



M 2 + a? -b 



+ — In 

4 



(M 2 + a 



2\ 2 



3tf 
4 



5 



(4tt) 



2 {6(6 - 4<?(F » + 2 [o(a - 4<7<y>i)) + pM] (M 2 + a 2 ) In [M 2 + a 2 ] 



+ 
+ 



a(4^ 1 )-a) + -(4^(F )-6)-pM 
a(4g( ( p l )-a)-^{Ag(F )-b)-pM 



(M 2 + a 2 + 6)E[M 2 + a 2 + 6] 
(M 2 + a 2 -6)hi[M 2 + a 2 -6]| > . (21) 



For the time being, what we have is a perturbative result for the effective potential. To 
actually obtain a nonperturbative result, we apply the optimization procedure. We had 
break the parameters My into ^-independent (ay) and ^-dependent (6y) parts and, with 
the help of (TT5|) . the optimized parameters turn out to be a i = a, 6n = 6, and p i2 = p. 
Upon application of the PMS, to find the optimized parameters a, 6 and p, we have to solve 
the three coupled equations 



av (1) 




dv {1) 




dv [l) 


da 


a=a,Q 


db 


b=b 


dp 







(25) 



P=P0 



at 5 = 1, and plug the optimized values ao, 6o and po into (I24jl . The following analytical 
solutions are encountered: 

a = 4a(<pi) = a 
6 = Ag{F ) = b , 

Po = = po . (26) 

The optimized parameters appear now as functions of the original coupling and fields, 
as expected. By plugging these results into (|24|) . all the S 1 terms vanish and the optimized 
potential can be cast under the form below: 



V. 



(i) 

eff 



57T 



16 5 2 <<pi) 2 ) 2 ln 



16g 2 (F o y 



+8^(F )(m 2 + 16^ 2 (<p 1 ) 2 )ln 



(m 2 + 16# 2 ((pi) 2 ) 2 
m 2 + 16g 2 (y 1 ) 2 + 4g(F ) 
m 2 + l6g 2 (<p 1 ) 2 -4g(F ) 



+16(7 2 (F ) 2 ln 



(m 2 + l6g 2 (<p 1 ) 2 ) 2 -lQg 2 (F )< 
p* 



48g 2 (F ) : 



(27) 
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The above result represents the Coleman- Weinberg-type potential 29|, |30) for the 



O'Raifeartaigh model 



13 



28} and it accounts for the sum of all one-loop diagrams. It 



is a nonperturbative result in that it takes into account all orders (actually infinite orders) 
in the original coupling constant. 



III. EFFECTIVE POTENTIAL IN THE LDE AT 0(5 2 ) 

We now present the order-5 2 results. At this order, we have one- and two-loop diagrams. 
Since the purpose of the present work is to be as clear as possible, the main part of this 
section is rather technical. All the results of this section are summarized in section IV. 



A. One Loop 

To analyze the one-loop diagrams we adopt the following strategy: there are 42 diagrams 
(plus some hermitian conjugates), and we separate them in distinct sets; each diagram 
belonging to a certain set has the same propagator structure in the loop and differs from 
the other diagrams of the set only by the external classical superfields or by insertions of a, 
b or p. Below, we show the expressions for these diagrams, set by set. We are going to use 
the notation Ji{6,6) for superspace integrals, which are listed in the Appendix. 

The first set is: 




; + h.c. 



FIG. 2: Diagrams g[ 1)S \ and 



(1)<5 2 



45V 



2^2 



d 4 kd% 2 

l2^F 
d 4 k 



$o(l)$o(2) 



1 = 



^(-fc)($i$i> 



+ h.c. 



(16) 2 " J (2tt) 4 
2<5V(^o) 2 



F{k)F{k)J 1 {d,8) + h.c. 



Kb 



(?7 2 ln?7 + — rfhirj — 2b) 



(28) 



-25 2 gb 



$O(1)02 



- A Dl{-k){^ x ) 



h.c. 



12 



(16) 2 

.Wo) 
Kb 



0^F(k)F(k)J 2 (6,6) + h.c. 



(r] 2 \nr} + - r] 2 \nr]- - 2b) . 



(29) 



:5 2 b 2 



d 4 kd 4 6 12 



4(16) 2 
S 2 b 2 , ^ 



(27T) 4 

16,4 / ^ F (k)Hk)MO,0) + h.c 



--DK-k)^) 



8Kb 

The second set is: 



(?7 2 ln?7 + — i] 2 \nrj — 26) 



+ h.c. 



(30) 



$1 






e 2 ■ e 2 { \ ; + /i.e. 



$]N /$ x 



FIG. 3: Diagrams and 



165V 



d 4 kd 4 9 12 

d 4 k 



$ (l)$ 1 (2) 



- iJ D 2 (fc)($ 1 $ 1 ) 



--Dli-k)^,) 



+ /i.e. 



/ (2^ FW " l^m(M)} +ft.a 



(16)= 

8Wo>W R+ _ Err) 



ft 



(31) 



-4S 2 ga 

_ 7 7TT \Q>\ b 

(16) 2 ' 1 
25 V 2 (F ) 

K 



d 4 kd 4 9 



12 



(2tt)< 



$o(l) 



--S?(*)<<M>i) 



--D 2 2 (-k)(%^) 



+ h.c. 



d 4 k 



^F(k) {A(k)J 6 (9,9) - \b\ 2 B(k)J 2 (6,6)} + h.c. 



(2ti 

(lnrj + — \nrj~) 



(32) 



?(1)5 2 



-45V 
45 2 g 



a\b\' 



dkdOu 2 

l2^ $l(1) ^ 2 
d 4 k 



\Dl(k)(^) 



~^(-*)<*i*o) 



+ /i.e. 



(16) 2 1 1 7 (2tt) 4 
2S 2 gab(ip 1 ) ,— . — _ 



F(fc) (A(£;) J" 7 (#, 0) - (^(^(M)} + ft.c. 



ft 



(lnr] + — Inr/ ) . 



(33) 
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g? 52 = 5 2 ab 
5 2 



d A kd 4 9 



(M 



12 el 



1 



1 



-_D 2 (-fc)($ $i> 



\a\ 2 \b\ 2 



(16)' 



d 4 k 



+ h.c. 

F(k) {A(k)J 9 (9,9) - \b\ 2 B(k)J 3 (9,9)} + h.c. 



The third set is: 





$1 



; + /i.e. 



FIG. 4: Diagrams , and g^* 2 . 



? (1)<5 2 



^*i(l)*i(2) 



8^ 



(27T) 

al 2 5 2 



-\p\{k)^ x ) 



~D 2 2 (-k)($ <S> ) 



+ h.c 



(16) 2 
k6 



(r] 2 ln?] + 



F(k)C(k)J 10 (9,9) + h.c. 
- ^EtT - 26) . 



/ d 4 kd 4 6\2 _ /1 , 



45 2 ^ 



I 1 212 

a a o 



(2 7 r) 4 



- 7J D 2 2 (-fc)($o$o) 



+ /i.e. 



(16) 2 1 1 7 (2tt) 4 
2S 2 ga 3 ( ( p 1 ) 



F(k)C(k)J 7 (9,9) + h.c. 



Kb 



[rj 2 \m] + — 7/ 2 ln?7 — 2b) 



1, 2 2 /• d A kd A 9 12 



-S z a 



(2tt)< 



--DKk)^) 



-±D 2 2 (-k)(%%) 



+ h.c. 



<r 2 



2 1 1 272 

a a o 



2(16) 2 
6 2 a A , 



Anb 



(27T) 4 

(r/ 2 Er/ + - ?7 2 Tn77~ - 2b) 



F{k)C(k)J 9 (9,9) + h.c. 
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The fourth set is: 




; + h.c. 



FIG. 5: Diagrams G^ , G^ and G^ . 



a 1 



(1)52 



85V _ 5 



d 4 kd% 2 
d 4 k 



$ 1 (1)$ 1 (2) 



+ /i.e. 



(16) 2 7 (2tt) 4 
+ |6| 4 J B(A;) J B(A:) l 7 12 (^^)} + / i .c. 
45Va 2 <^> 2 



{A(A;)A(fc) Jn(0, 0) - 2|6| 2 A(A;)B(fc) t 7 1 o(0, 0) 



[461nr/ 2 - V + 2&)lnr/ + + (r? 2 - 2&)ln7/~ + 2b] . 



(38) 



»12 



i (i kd 4 9\o , . 
-45V / ^ 



45 2 # 



a| 2 a 



(2^) 4 
d 4 /c 



--D 2 2 (-k)(^ ) 



+ h.c. 



(16) 2 ' 1 J (2vr) 4 
+ |6| 4 B(A;)5(A;)J 8 (M)} + /i.c. 
25 V 3 (y^) [AlT - 



{A(k)A(k)J 13 (9, 9) - 2\b\ 2 A{k)B{k)J 7 {9, 9) 



nb 



[4b\nr] 2 - {if + 26)lnr/ + + (r? 2 - 2b)\mf + 26] 



(39) 



« 52 = ~/a 2 



ctkc£9u 
d 4 k 



--Dl(k)(^ ) 



^D 2 2 (-k)(^ ) 



h.c. 



{A(k)A(k)J u (9, 9) - 2\b\ 2 A(k)B(k)J 9 (9, 



2(16) 2 1 1 J (2tt) 4 
+ |6| 4 J B(A;) J B(A:)J3(^^)} + /i.e. 
5 2 o 4 

[4&ln7/ 2 - (?? 2 + 26)ln7/ + + (r? 2 - 2&)lnV + 2b] . 



(40) 



The fifth set is: 





9 2 + /i.e. ; fl 2 / M 2 

<i>i — y|> 1 



FIG. 6: Diagrams G { £ & \ G^ and . 
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;(1)5 2 
'14 



85V 
8S 2 g 



-±5?(*)<<Mi> 



-±Z> 2 2 (-fc)<*i<l>i> 



16 

(16)= 



■B(k)B(k)J 17 (6,6) 



k ye 
2r) 2 n 2 (F ) 2 

+ y : 0/ [4felnr/ 2 - (?7 2 + 26)lnr? + + (r/ 2 - 2&)lnrf + 26] . 



k6 



j(1)5 2 
'l5 



+ /i.c 



25*9 
16 



5 /^{^(^m8(^^) + 2^Wi3(A:)J 2 (^^) 



■B(k)B(k)J 19 (9, 9) \+h.c. 



(16) : 

4^W)/l_- 2 
k \e 



k6 



[46I1177 2 - (r? 2 + 26)ln7]+ + (r? 2 - 2&)ln?r + 2b] 



^16 



1 r*kd%> 

2 11 J (2tt) 4 
5 2 



-±£>2(-*)<*i*i> 



2(16) 



|6| 



d A k 



B(k)B(k)J 21 (9,9) 



E(k)E(k)J 20 (9, 9) + 2^E(k)B(k)J 3 (9, 9) 



2 ( 2 16)2 

— \nr] 2 



2k \ e 
<5 2 6 2 



+ 



8k& 



[461nr] 2 - {rf + 26)lnr] + + (r/ 2 - 2b)lnrj~ + 2b] 



The sixth set is: 



*1 




$1 N 



$ i . $1 

$1 $ 




*1 



$1 



$1 




$1 $1/— >v$0 



2 ; 9 2 



U; + h.c 



FIG. 7: Diagrams Q^ S \ Q { ^ S \ and g$* 2 . 
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G {1)s2 



165V / d -^Mi)H2) 



■ah 



- i £^(-fc)($ $i> 



(2tt)4 

^C(fc) <j E(k)J 22 (0, 9) + J^.fl(fc) J 23 (0, 0) |> + /i.e. 



+ /i.e. 



(16)2 

8*ya<Fo)M /t — + i — 



ft 



(lnr? 4 " — In?] ) 



(44) 



»18 



_ / d 4 kd 4 9\ 2 / 1N 
-4*V / -(^ $ o(l) 



- i ^(fc)($ 1 $ 1 > 



--Z>a(-fc)<So*i> 



+ /i.e. 



45 V ,27 

'a 2 & 



(16) 2 

_ 9 — (__o)_ ^j n ^+ _ 



{ . H)A C{k) <j E(k)J 6 (9,9) + ^5(^(0,0) >> + h.c. 



(45) 



/-(i)5 2 Art i f d A kd A 9i 2 - 
^19 = "^V y (2tt) 4 l( ^ 



-£> 2 (£;)(<ii$i 



1 



D 2 2 (-fc)($!$o> 



+ h.c. 



45 2 g 



a\b\' 



d 4 k 



(16) 2 1 1 J (27T) 
25 2 gab(<p 1 ) 



4 C(fc) <J E(k)J 25 (9, 9) + -^B{k)J 2& {9 } 9) \ + /i.e. 



16 



ft 



(46) 



^20 



<5 2 a6 7 d ^ 6,12 fl2 



(2tt)< 



|a| 2 |6| 2 



(i 4 A; 



(16)2' 1 1 1 J (2n 
5 2 a 2 b ^ + ^ 



-\n\{-k)^) 

|6| 5 



+ /i.e. 



4 C(fc) <j E(k)J 9 (9, 9) + -^B(k)J 27 (9, 9) \ + /i.e. 



16 



2/t 

The seventh set is: 



(ln?7 + — ln?7 ) 



(47) 



$1 




$1 $i 



$1 




$ $ /— "X^o 
+ h.c. ; 

$1 <!>! 



FIG. 8: Diagrams , and . 



»21 



16^V / ^^(1)^(2) 



-^(*)<*i*i> 



^ 2 2 (-fc)(l»o$o) 



17 



16 

+(k 2 + \M\ 2 )^B{k)A{k) J 30 (9, 6) + 



d 4 k 

I' 

16 



{(A: 2 + \M\ 2 )E(k)A{k)J 28 (9,9) + \a\ 2 \b\ 2 E(k)B(k) J 29 (9, 9) 

a\ 2 \b\ 4 



16 



-B(k)B(k)J 31 (9,9) 



{4b(27] 2 - M 2 )\nrf + \a 2 if + 2r] + (M 2 - r} + )]\nr] 



+ 



-[a 2 r] 2 + 2rf{M 2 - rf )\\mf - 2b(a 2 + 2M 2 - 2r] 2 )} 



(48) 



7 (1)<5 2 
'22 



-45 2 #a 
45 V 



d 4 kd 4 9 12 
d 4 k 



-±d?(ao<<Mi) 



--£ 2 (-£;)<i> cJ> ) 



+ /i.e. 



(2tt)' 

!&! 

16 



{(A; 2 + |M| 2 )£(fc)A(AOJi 3 (M) + |a| 2 |6| 2 £(£;)5(fc)J32(M) 
'a| 2 |6| 4 



+ (A: 2 + \M\ 2 )^B(k)A(k) J 7 (9, 9) + 



16 



-B(k)B(k)J 8 (9,9) } + h.c. 



25 2 ga((p 1 ) 
Kb 



{Ab{2r] 2 - M 2 )\nrf + [a 2 r/ 2 + 2r] + (M 2 - r/ + )]hV 



'2.3 



^ = 5 2 |a| 2 
lay 



-[aV + 2r]-(M 2 - rf )]lnrT - 2b(a 2 + 2M 2 - 2r] 2 )} 
d 4 kd 4 9 12 ~ 



(49) 



(2n) 4 
d 4 k 



\D\{k){^i) 



- A D 2 2 {-k){%%) 



{(k 2 + \M\ 2 )E(k)A(k)J 33 (9,9) + \a\ 2 \b\ 2 E(k)B(k)J u (6, 9) 



16' 

+(k 2 + \M\ 2 )^B{k)A{k) J 35 (9,9) + 



a| 2 |6 | 4 
16 



16 



-B(k)B(k)J 36 (9,9) 



5 2 a 2 
4Kb 



{4b(2r] 2 - M 2 )\nrf + [a 2 rf + 2r] + (M 2 - 7?+)]hV 



-[a 2 ?? 2 + 2r]~(M 2 - r]-)]\nr]- - 2b(a 2 + 2M 2 - 2rf )} . 
The eighth set is: 

^ox—n^i ^ ^ 4> /— $0/— 



(50) 




4> 




+ h.c. 



4V 



FIG. 9: Diagrams , and ^e^- 
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165V 
165V 



d 4 kd 4 e 12 



*i(l)$i(2) 



■±D?(fcX*i*o> 



±D|(-*X*i*o> 



\a\ 2 \b\ 2 



d 4 k 



(16) 3 1 1 1 1 J (2tt) 4 



C(k)C(k)J 37 (9,9) 



Kb 



(rj 2 lnrj + - r/ 2 ln7/ _ - 26) . 
18 



(51) 



»25 



= -4(5 ' ga 
45 2 g 



d A kd A 6 12 



-^(*)(^o) 



1 



+ h 



a\ 2 d\b\ 2 



d A k 



(16) 3 ' 1 1 1 J (2n) A 
25 2 ga 3 (<p 1 ) . 2 - 



c(k)c(k)j 38 (e,e) + h. c . 



nb 



{f?\mf - r] 2 \nr]~ - 2b) 



'26 



= 51a 



2 I „ 1 2 



d 4 kd A 12 



--5?(fc)<*i*o> 



-Z>2(-fc)<*i*o) 



|4|7 |2 



(i 4 A; 
(2^ 



(16)3' 1 1 1 

5 2 a 4 — — 

(i] 2 \nrj + — r) 2 \nr] — 2b) . 



C(k)C(k)J 39 (9,9) 



The ninth set is: 




$ 2 



32 / j ; -i /*.<•. 



FIG. 10: Diagrams g^ 5 ' and . 



»27 



i d 4 kd 4 6\2 , , 

- 45 V / /O $ (>(1) 



(2tt) 4 
d 4 £; 



-^(fc)($i$ 2 > 



(16)^ M6 7 (2rr) A 



F(fc) {4(fc) J 6 (0, 0) - \b\ 2 B(k)J 2 (9, §)} + /i.e. 



25 2 g P M(F ) - 

K ' 



^28 



5 2 bp 



«12 



0? 



- i 5 2 (^)($ 1 $ 2 ) 



d 4 A; 



(16) 2r 11 J (2tt) 4 
^ (Er^-Err) 



F(k) {A(k)J 9 (9, i 



--DK-fc)^^) +/I.C. 

•)-|6| 2 S(fc) t 7 3 (0,0)} + /i.c 
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The tenth set is: 



$1 




$1 



; + h.c. 



FIG. 11: Diagrams G29" and Q 30 



$ 2 



M)S 2 

'29 



-46 2 gp 
4S 2 g 



d 4 kd% 2 



~^(fc)($l$ 2 > 



-^(-*)<*i*o) 



(2vr) 4 

/ 7 ^l{^( A; )^)^i3(^^-2|6| 2 A(A:) J B(A;)J 7 (^^ 



(16) 2 ^ 7 (2tt) 4 
+ \b\ 4 B(k)B(k)J s (6,9)} + h.c 
25 2 gpMa(tp 1 ) 



Kb 



[4&lnr/ 2 - (rf + 2b)\nr] + + (rf - 2b)\nrj~ + 2b] . 



»30 



£ 2 ap 



d 4 kd 4 9 12 



--£> 2 (£;)<$i<t>2> 



-l5|(-fc)($i$ > 



+ /i.e. 



P M|a| 2 / 7 ^ I {A(A:)A(A;)Ji 4 (^^)-2|b| 2 A(A;) J B(A:)J 9 (^^) 



(16) 2r 1 1 J (2vr) 4 
+|fc| 4 £(fc)£(£;)J36(M)} + /*-c. 
5 2 pMa 2 r - 



2nb 

The eleventh set is: 



[4&I1177 2 - (r? 2 + 26)lnr] + + (rj 2 - 2b)\mj~ + 2b] 




$2 $0/— 

x x ; + /i.e. 

>iV_y$ 1 



$1 



FIG. 12: Diagrams 0^ and 



M)S 2 
'31 



d 4 £;d 4 # 



12 



(2tt) 4 



-^(fc)<<Ml> 



(16) 2K 7 (2tt) 4 
25 2 gpMa{ip 1 ) 
Kb 



F(k)C(k)J 7 (9,9) + h.c. 
(?7 2 ln?7 + — ?7 2 lnr]~ — 26) . 



•^ 2 (-fc)($ 2 $ ) 
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7 (1)<5 2 
'32 



5 ap 



(16) s 



d 4 kd 4 e 12 

pM|a| 2 6 2 



-lD 2 2 (-fc)($ 2 $ ) 



+ /i.e. 



c/ 4 A; 
(2^ 



F(k)C{k)J 9 {9,9) + h.c. 



$ 2 pMa 2 , 2 — + 2 — _ , 
2 & lnr ? - ^/ Inr; - 26) 



The twelfth set is: 



civ 




32 




$1 



; + h.c. 
FIG. 13: Diagrams and 



c (1)<52 

»33 



»34 



p^o(l) 



+ h 



16 



d^k 



F(k) t E(k)J 18 (9, 9) + U-B(k)J 2 (6, 9) \ + h.c. 



25 2 gpM(F ) ,- . - 



K 



(lnr? 4 " — Inrj ) 



5 2 bp 



«12 



~pM\b\ 2 



-el 
d 4 k 



\Dl{k){^ 2 ) 



F{k) E(k)J, 



20\ 



S 2 pMb 
2k 



(ln^ + — Inr] ) 



16 



16 



361 



+ /l.C. 

+ /l.C 



The thirteenth set is: 




<J>i^ — ^(1)2 $1^- / <I>2 



; + h.c. 



FIG. 14: Diagrams e^" 5 ' and . 



M)S 2 

'35 



/ 77^1 $ l(!) 



(2tt)< 



■±5?(fc)<*l*2> 



±Z?2(-*)<*1*0> 



+ /l 



d 4 A; 



(16) 2r 1 1 J (2^) 4 
25 2 gpMa(y 1 ) . 2 - 



F(k)C(k)J 7 (9,9) + h.c. 



Kb 



(rj 2 lnrj + - rj 2 \nr]~ - 2b) . 
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M)S 2 
'36 



S 2 ap 
5 2 



d A kd 4 e 12 



r P M\a\ 2 \b\ 2 



-\D\{k){^ 2 ) 



d A k 



--I?2(-*)<*i*o> 



+ h.c. 



(:IG) --- / { ,^F(k)c(k)j 9 (e,e) + h. c . 



2 & I 7 ? ln V ~ V In?? - 26) 



The fourteenth set is: 



(63) 




$2 $0X~\$2 
! $]V /$ 1 



+ /l.C. 



FIG. 15: Diagrams 0§ )tf2 and g$ fi . 



r d A kd A 9 

J 



45^ 
16 



pMd 



^(l) 
d A k 



--5?(fc)<*i*i 



-^ 2 2 (-fc)($2<t>o) 



+ /i.e. 



(2tt)< 



{-£(£;),4(£;)j4o(M) + \b\ 2 E(k)B(k)J 32 (9,9) 



-\^-B(k)A(k)J 7 (6, 9) + J" 8 (0, 0) |> + /i.e. 

25 2 gpMa{(pi) 



Kb 



[Ab\mf - (?? 2 + 2b)\m] + + (r? 2 - 2b)\n.-q~ + 26] 



(64) 



^38 = 0/9 



d A kd A 9 12 
~J2^ 



±D 2 2 (-k)(<j> 2 %) 



h.c. 



^pM\a\ 2 J -0^ {-E(k)A(k)J 33 (9, 9) + \b\ 2 E(k)B(k)J u (9, 9) 

-\^B(k)A(k)J 35 (9, 9) + \^B(k)B(k)J 36 (9, 9)} + h.c. 
5 2 pMa 2 



2Kb 



[Ab\nr] 2 - (rf + 2b)\m] + + (rj 2 - 2b)\nrj- + 26] 



(65) 



The last four diagrams have unique propagator structure in the loop and are not divided 
into sets. 

The thirtieth ninth diagram is: 



G 0)P 

»39 



d A kd A 9 



12 



(2tt)< 



- A Dl{k){^ 2 ) 



\D 2 2 {-k) (Ml) 



+ /l.C. 



A 2 f d 4 k 

6 -p 2 \M\ 2 b 2 ' : " " 



(2tt)^ 



2(16) 2 

= ^(^W-^T-26) . 



C(k)F(k)J 9 (9,9) + h.c. 



(66) 
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The fortieth diagram is: 



$ip- 
X x + /i.e. 



-/<£> 2 

FIG. 16: Diagram Q m 



(1)5 2 



•p., 

FIG. 17: Diagram . 



»40 



5 2 |P| 2 



d 4 A;(2 4 ( 



12 



(2tt)< 



—^(-A^iSi) 



r) 2 r rl 4 k ( - \b\ 2 

= Yq \p? J |(fc 2 + |a| 2 )^(A:)£;(A:)j733(^, ^) + + |a| 2 ) ^-^(A:)i?(A:) ^C^, ^) 



|M| 2 |fr | 4 
16 



+ \M\ 2 \b\ 2 B(k)E(k) J 34 (9,9) + 



-B(k)B(k)J 36 (9,9) 



S 2 p 2 
4Kb 



{Ab(2r] 2 - a 2 )\nr] 2 - [27] + (r] + - a 2 ) - MV]lmf 



+ [2r/-(7y _ - a 2 ) - M 2 rf]\nr]- + 26(2?? 2 - 2a 2 - M 2 )} . 



(67) 



The fortieth first diagram is: 



$ 2 V — /$ 2 



+ /i.e. 



FIG. 18: Diagram ^J^. 



f 41 



1 2 2 /" d 4 kd 4 9 l2 
-dp ' 



(2tt)< 



-±D*{k)(Q 2 *i) 



+ h.c. 



; P 2 M 2 



d A k 



2(16) 2r J (2tt) 4 
+ \b\ 4 B(k)B(k)J 3 (9,9)} + h.c. 
S 2 p 2 M 2 r - 



{A(k)A(k)J 14 (e, 9) - 2\b\ 2 A(k)B(k)J 9 (8, 9) 



4Kb 



[4b\nr] 2 - (if + 2b)\m] + + [rf - 2b)\mf + 2b] 



(68) 
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<J> 2 V /<j> 



FIG. 19: Diagram . 



Finally, the fortieth second diagram is: 



t?42 = IPl 



d 4 kd 4 9 12 



^|plW|6| a /^(*W)^ 

— ^— - — (r) 2 \nri + — r] 2 \nr)~ — 2b) . 
Ako 



20 1 



(69) 



The diagrams Q[^ 5 , Q[^ S and Q[q° of the fifth set give divergent contributions and we 
will discuss their renormalization later. 



,(1)5 2 



B. Two Loops 



The two-loop diagrams are shown in Fig. [2UJ 

+ h.c. + + h.c. + 



$ $1 




$ $1 




$1 $c 



+ 



$1 <E>1 $1 $1 




FIG. 20: Diagrams Q^ 5 \ Q^ S \ and gf & \ 



In order to calculated the two-loop diagrams we are going to use technique developed in 



31] . The contribution of the first diagram is given by 



(2)<5 2 



. c n o f d pd kd 9 \2 

4:0 q / ; r-T 

J W 

S 2 Q 2 f d 4 pd A k 



-~Dl(k) 



(70) 



(16) 3 7 (2tt) 5 

with g = (k — p) and the integral I\(9,9) is given in Appendix A. Plugging X\{6,6) into 
7D|) . we obtain: 



3} 



(2)5 2 



-45Va 2 £> 3 



d 4 pd 4 k 
(2vr) 8 



C(p)F(k)F(q)p 2 



(71) 
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To handle with the integral above and the other momentum space integrals that will ap- 
pear in the following, we have adopted the strategy: for each of them, we split the integrand 
with the help of the method of partial fraction decomposition and write each integral as the 
sum of other integrals with just three terms in the denominator. The remaining integrals 
are well known in the literature, and we use the results of the references 32|-l34j| to compute 
them. 

From now on, we define rf = m 2 + a 2 , rf 1 = m 2 + a 2 ±b and adopt the same notation of 



references 



32 



34| for the integrals I(x, y, z) , J(x, y) and J(x): 
kJ(x) = h x (inx — l) , 



k 2 J(x, y) = xy 



— H — (2 — \nx — In y) + (l — lnx — \ny + lnxlny) 



(72) 
(73) 



n 2 I(x,y,z) 



1 - L^j - i {L 2 - 6L1 + (y+z-x)lnyhiz 



c 

e V 2 "7 2 

+ (z+x — y)\nz\nx + (y+x — z)lny\nx + £(x, y, z) + c [7 + C(2)]} , (74) 



where 



k = (An) 2 , 



c 

\nX 

i(?,y,z) 
s 

Li 2 (z) 



x + y + z , 



In 



X 

7 2 



+ 7 — In 4-7T 



xm x + ym y + zln z , 
S 



2 In Z + X - V - - In - + V - X - - - In - In - 



-2Li, 



2z 

z + x — y — S 
2z 
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-2Li 2 



z z 
z + y — x — S 
2z 



\/ x 2 + y 2 + z 2 — 2xy — 2yz — 2zx 
_ r ln(l-t) 

Jo t 



-dt (dilogarithm function) . 



7T 



So, the total contribution of the first two-loop diagram is: 
g[ 2)5 = — - — [l(r] + ,rj + ,r] + )-3I(r] + ,ri + ,r]~) + 3I(r] + ,r]~,ri')-I(r]',ri'. i r]')] 
The contribution of the second diagram is given by 



£ (2)<J2 = 85 2 g 2 [ d *P d4kd4dl2 



(27T) 



^D 2 2 (k)(^i) 



(75) 



^Dl{q)Dl{-q){^ x ) 
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OA 2 /! 2 r rfind^k 

JW a2b J W~ {^WW^M - 2b 2 A(p)B(k)F(q)l 3 (9,< 
+b 4 B(p)B(k)F(q)Z 4 (9j)} . 



(76) 



Plugging Z 2 



and X 4 ( 



into (1761), we obtain: 



,(2)5 2 



-85 2 g 2 a 2 b 3 { 2 



d ' pdAk A(p)B(k)F(q)p 2 + b 2 fd4pd ' k 



(2tt) 8 v ' v ^ J (2tt) 8 

Decomposing the momentum space integrals we get the contribution: 



B{p)B{k)F{q) . (77) 



5 2 g 2 a 2 [-4/(?] 2 , r? 2 , 77+) + 4I(r) 2 , rj 2 , ry") + 1(77+ 77+ 77" 
+/(?7 + , ?7 + , r]~) - I(rj + , r}~,r]~) - I(?T, 77", 77" 



(78) 



The contribution of the third diagram is given by 



irX 2 2 I'drpd 4 kd 4 9 12 
16o g 



~D 2 (k){^i) 



(2tt) 8 

., rJ , ^ pd k [ , . -. I o 

a b ■ 



16 



D?(g)Df(-g)<$ ^i> 



5 V 2,2 /" ^V 4 ^ 



8 C(p)C(q) { E(k)l 5 (9,9) + -b 2 B(k)l 6 (9, 



(16) 3 J (2n) 8 _Vf/ " w ^-^--v-'-' 1 i 6 
Plugging Z 5 (0, 0) and Zg^, #) into (f79"|) . we obtain: 

Decomposing the momentum space integrals we have: 



(79) 



C(p)B(k)C(q)p 2 q 2 



(80) 



,(2)«5 2 



2S 2 g 2 a 2 [l(t] + ,r] + ,ri + ) - I(r] + ,r] + ,rj )-I(r] + ,r] ,77 ) + I(r} ,77 ,77 )] . (81) 



The contribution of the fourth diagram is given by 



MVS 2 8i2 2/ drptrktr On 
Gi =86g J (2vr)B 

d 4 pd A k 
(2tt) 8 



4 



16 



r) 2 ( g )D 2 2 (-g)($ $o> 



85 V , 
(16) 2 
1 



B(q) E(p)E(k)l 7 (0,6) + -b 2 E(p)B(k)l 8 {9,> 



+ 



; b 4 B(p)B(k)U0, 



(16)2 

SS 2 g 2 fd A pd A k 



my2 / rH)8 -A(q)(q 2 + m 2 )^E(p)E(k)l 10 ( 



\b 2 E{p)B{k)lii{9,9) 



(16) 



; b 4 B(p)B(k)l 12 (9,9) 



(82) 
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Plugging X 7 (9, 6) - T 12 (6, 9) into §2]), we get: 



,(2)<5 2 



85 2 g 2 b 2 { a 2 



d 4 pd 4 k 



E(p)E(k)B(q) + 2a 2 b 2 



d 4 pd 4 k 



E{p)B{k)B{q) 



+a 2 b 4 



d 4 pd 4 k 
(2tt) 8 



B(p)B(k)B(q)-2 



d 4 pd 4 k 



E(p)B(k)A(q)k 2 q 2 



-b 2 i d ^B{p)B{k)A{ q )q 4 -2m 2 



-m b 



(27T) 
2,2 /" 



(2tt) 



E(p)B(k)A(q)k 2 



(2nf 



B{p)B{k)A{q)q 2 



(83) 



Decomposing the momentum space integrals, we obtain the total contribution of the fourth 
diagram: 



,(2)5 2 



5 2 g 2 { 8a 



-i(v 2 , v 2 , v 2 ) + Xhv 2 , v 2 , v + ) - Xhv 2 , v 2 , V ) 



+8m 



-2I( V 2 , rj 2 , 0) + ^/(t? 2 , rj + , 0) + \ r,', 0) 



+a 2 [I(r] + , f] + , r] + ) + 3/(?7 + , rf) + 3i"(r? + , r)~,r}~) + I(t]~,r)~, 7/~)] 
+2 [-4J(?7 2 , ?7 2 ) + 4J(?7 2 , 77+) + 4J(?? 2 , tj") 
-J( V + , V + )-2J( V + , v -)-J(rj-, rj')]} ■ 

(2)<5 2 



(84) 



Unlike the first three two-loop diagrams, Q\ gives a divergent contribution, however, 
its renormalization is trivial, since it is a vacuum diagram. 



C. Regularization and renormalization 

The divergent diagrams of order 5 2 are: 
85V A r- 



»14 



»15 



sir 2 



k \ e 

2„2 



lnr^ 



2Pg 
Kb 

2S 2 gb fl 



[4b\n v 2 -(r] 2 + 2b)\nr} + + (rj 2 - 2&)lmT+2&] j ^ d 4 #$ $o ; 



In?/ 



[4b\m 1 2 -(r ] 2 + 2b)\nT] + + (r] 2 -2b)\nT]- + 2b\ ^Jd 4 66 2 % + h.c. 



2Kb 

S 2 b 2 I 1 



2k \e 



hurj' 



i5) 



(86) 



+— [4&E77 2 - (?? 2 + 2b)Erj + + (?? 2 - 26)E77" + 2b ]\ / dWfl 2 • (87) 
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The diagram Q^ s is a vacuum diagram and thus its renormalization is trivial. To 

2 

w< 

85 2 g 2 



renormalize the divergent term in Q[^ S we introduce the counterterm 



K€ 



d A e$ R®0R , (88) 



and, for the divergent term in 0[^ 5 , we introduce the counterterm 

^ / d A ee 2 % R + h.c. = ^ / d 2 e% R + h.c. (89) 

Ke J ne J—. 

As mentioned in the previous section and showed in 14J, plugging the solution for the op- 
timized parameter b in ( 189]) . only the Kahler potential is actually renormalized, in agreement 
with the nonrenormalization theorem. 



IV. SUMMARY OF THE RESULTS AND NUMERICAL ANALYSIS 



In this last section we summarize the perturbative results for the effective potential up 
to the order 5 2 . We also derive nonperturbative corrections to the effective potential by 
implementing the PMS criterion numerically. Up to this order, the renormalized effective 
potential can be written as 



' // - V ef7 + V eff + V eff + V e// • 



Below, we write separately the terms of the effective potential V e // as follows: 
The tree level potential is given by: 

with £ < 0. 

The vacuum diagram of order 5° is: 
,s° _ 



V. 



(4{i (M2+ ° 2)2ln 



(M 2 + a 2 



+ U -{M 2 + a 2 )\n 



M 2 + a 2 + b 



+ — In 

4 



(M 2 + a 2 f-b 



M 2 + a 2 -b I ^ 

The one- loop contribution for the effective potential at 0(S l ) is given by: 
5 



3^ 



v 5 ff 



(An) 2 
+ 



{b(b - Ag(F }) + 2 [a(a - Agfa)) + pM](M 2 + a 2 ) In [M 2 + a 2 ] 



1 



+ 



a(4g(ip 1 )-a)+-(4g(F )-b)-pM 
a(Ag{ Vl ) - a )-~(4g{F )-b) -pM 



(M 2 + a 2 + b)hi[M 2 + a 2 + b} 
(M 2 + a 2 -b) hi [M 2 + a 2 -b] 



(90) 



(91) 



(92) 



(93) 
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At the 0(S 2 ), we separate the contributions of one- and two-loop diagrams: 



V S 2 _ y<5 2 , y<5 2 

y eff ~ V eff(I) + V eff(II) , 



(94) 



where 



V eff(I) 



a - Agfa}) 2 (M 2 + 3a 2 ) + apM(a - Agfa)) + p 2 (3M 2 + a 2 )] 



(An) 2 
xE (M 2 + a 2 ) 

- I \(a - Agfa)) [(a - Agfa})(M 2 + 3a 2 + b) + 2a(b - Ag(F }) + AapM] 



In (M 2 + a 2 + b) 



+ \{b- 4<7(F » [{b - 4 (F » + ApM] + p 2 (3M 2 + a 2 + 6) 
- 1 [(a - Agfa)) [(a - Agfa))(M 2 + 3a 2 - b) - 2a(b - Ag(F )) + AapM] 
+ \{b- Ag(F )) [(& - Ag(F )) - ApM] + p 2 (3M 2 + a 2 - b) 



x In (M 2 + a 2 - b) } , 



(95) 



and 



V, 



5 2 



eff (II) 



{ (M 2 + a 2 ) ^ 4(M2 + ^ + 6(M ' + ^ " b)(M2 ~ 2 ° 2) ] E2(M2 + fl2) 

- ^ll^X^ K M ' + « 2 )( 6M2 + 17 « 2 ) + 2& (M 2 + 3a 2 )] hT 2 (M 2 + a 2 + 6) 
(M + er) 

-3(M 2 + a 2 - b)(2M 2 + 3a 2 - 26)E 2 (M 2 + a 2 - b) 
+8(M 2 + a 2 + o)(2a 2 - 6)ln"(M 2 + a 2 )IK(M 2 

n 9 T \ 



+ a 2 + b) 



8a 2 b(M 2 + a 2 - b)— , - 2 , r , lf2 2 ,. 
f ( 1 m2 + q2) J ln(M 2 + a 2 )ln(M 2 + a 2 - 6) 

-2(M 2 + a 2 + 6)(2M 2 + 3o 2 - 26)E(M 2 + a 2 + 6)E(M 2 + a 2 - 6) 

8M 2 6 L ( (M 2 + a 2 ) 2 \ . „ 2 . n (M 2 + a 2 + b^ 

bln TaF? oJ I? - M 2 + a 2 In — — - 

\(M 2 + a 2 ) 2 -b 2 J \M 2 + a 2 -b 



ln& 



(M 2 + a 2 )^ 
-24(M 2 + a 2 )(2M 2 + 3a 2 )hT(M 2 + a 2 ) 

+ 12(M 2 + a 2 + 6)(2M 2 + 3a 2 + 2o)E(M 2 + a 2 + 6) 

+ 12(M 2 + a 2 - 6)(2M 2 + 3a 2 - 26)E(M 2 + a 2 - 6) 

+a 2 [4£(M 2 + a 2 , M 2 + o 2 , M 2 + a 2 )-3£(M 2 + a 2 + b, M 2 + a 2 + b, M 2 + a 2 + b) 
-C(M 2 + a 2 + b, M 2 + a 2 -b, M 2 + a 2 -b)] 

+ 4a2 % + + a ^ b) [e(M 2 + a 2 ,M 2 + a 2 ,M 2 + a 2 + 6) 
- £(M 2 + a 2 ,M 2 + a 2 ,M 2 + a 2 -b)] 
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8M 2 b 



(M 2 + 



(M 2 + a 2 + b)U 2 



' M 2 + a 2 
M 2 + a 2 + b 



+ {M 2 + a 2 - &)Li : 



M 2 + a 2 -b 
M 2 + a 2 



5Qb 2 --bn 2 (2M 2 + b) 



(96) 



A. Optimization Procedure 

Let us now describe the numerical results obtained by implementing the PMS criterion. 
Our numerical approach consists in solving a system of coupled equations, to determine the 
solutions for the optimal parameters a, b and p that satisfy the PMS criterion and minimize 
the effective potential up to the order 5 2 . We start comparing our numerical results with 
the analytical results derived at order 5 1 and implement the optimization at order S 2 . 

The numerical implementation was performed using Mathematica 35]. When we apply 
the PMS criterion to the effective potential there are three coupled equations to solve: 

dVeff 



da 
db 
dp 



, 
, 
. 



(97) 



We use the effective potential evaluated up to the order 5 2 , which is given by Eq. fl90|) . 
The criterion that we have stablished to choose the optimum roots is to follow the same 
family of roots and work with the roots that minimize the effective potential. Although it 
is not possible to derive analytical solutions to the PMS equations up to the order 5 2 , it can 
be seen that the solutions can be written as 



a = AgM + hAig,®,®) , 
b = 4g(F } + hB(g,<S>,<S>) , 
p = + HC(g,$,$) , 



(98) 



where A, B and C are corrections to the order 5 1 solutions. This naturally brings nonlinear 
g contributions and generates nonperturbative results that go beyond the one- and two-loop 



results derived in 



13 



3l|. 



The parameters were chosen so that spontaneous supersymmetry breaking appears in the 
O'Raifeartaigh model. We choose £ = —10 and m = 10 and perform one rescaling in all 
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quantities in terms of the renormalization scale \i or, in other words, our quantities are given 
in unities of p. In Fig. [21] we show the tree level effective potential for the O'Raifeartaigh 
model, and we can see that SUSY is spontaneously broken, since the value of the effective 
potential at the minimum is different from zero. 
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FIG. 21: The results for tree level effective potential for different values of the coupling constant 
g. Parameters: £ = —10.0, m = 10.0. 



In a previous work 13] was shown in detail that different optimization procedures FAC 
(Fastest Apparent Convergence) 2 and PMS give the same result for the optimal parameters 
a, b and p at 0(5 l ), and at this order was possible to implement these two optimization 
procedures analytically. 

The effective potential obtained in this case is shown in Fig. (22J where we show the results 
obtained up to the order 5 l using the PMS and FAC criteria. Our analytical and numerical 
results including only Old 1 ) contributions to the effective potential have shown that the 
solution that minimizes the effective potential is p = 0. In this figure we apply the PMS 
criterion numerically at Old 1 ) and compare with the effective potential analytically evaluated 
with the FAC and PMS criteria. We change the value of the coupling constant g and we can 
see that for different optimization procedures we obtain the same effective potential, as it 
should be, based on the analytical results obtained by the different optimization procedures. 

We start with the comparison between the nonperturbative effective potential up to the 
and the 0(5 2 ), allowing us to gauge the performance of each optimization procedure, 
regarding both reliability and indications of the convergence of the method. 



2 for a description of the optimization procedure FAC please see Ref. [3] 
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g = 0.10 g = 0.50 




-4 -2 2 4 -4 -2 2 4 

tfil <Pl 



FIG. 22: The results for V e ff(ipi) using the optimization procedures FAC and PMS up to the 0(5 l ) 
in LDE for different values of the coupling constant g. Parameters: £ = —10.0, m = 10.0. 





V cS fa = 0] 




PMS 0{5 l ) 


FAC 0(5 X ) 


PMS 0(5 2 ) 


g = 0.01 


100.002 


100.002 


100.002 


5 = 0.10 


100.233 


100.233 


100.231 


g = 0.50 


105.828 


105.828 


104.544 


g = 0.75 


113.101 


113.101 


107.394 


g = 1.00 


123.260 


123.260 


107.865 


g = 1.20 


133.448 


133.448 


105.631 



TABLE I: Dependence of the minimum of effective potential with the coupling constant g including 
the corrections at 0(5 2 ). In this table we show the analytical results using FAC and numerical 
results for PMS. Parameters are the same of Fig. [22] 

As a result, we find that at the 0(S 2 ) the solution that minimizes the effective potential 
is again p = 0. We can see in Table |T] the results for the value of the minimum of the 
effective potential, and we note that the results at the 0(5 2 ) are very similar to the results 
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at the O^ 1 ) up to g ~ 0.5, indicating that for the parameters used here, g ~ 0.5 can be 
a regime of strong coupling. To our knowledge this is the first work in the literature that 
performs the evaluation beyond the one-loop approximation of the effective potential for the 
O'Raifeartaigh model (with spontaneous SUSY breaking). For the parameters considered 
in this study we can see strong indications that the nonperturbative method of LDE is 
appropriate to deal with the O'Raifeartaigh model. The results at the 0{5 1 ) are the same 
as those at the 0(5 2 ) up to some value of the coupling constant, indicating some convergence 
of the nonperturbative method. 

V. CONCLUDING REMARKS 

The investigation we have pursued in this paper and the results of our explicit supergraph 
evaluations confirm that superspace techniques, even if SUSY is spontaneously or explicitly 
broken, can be consistently combined with the LDE approach to compute higher order cor- 
rections to effective potentials in the framework of supersymmetric models. We just highlight 
that we are here bound to the minimal O'Raifeartaigh model, but the extension to general- 
ized (Wess-Zumino type) matter models can naturally be pushed forward. Fayet-Iliopoulos 
D-term breaking models are also very interesting to be reassessed with the approach we have 
adopted in the present work. It is true that these first attempts to go through higher orders 
with the LDE procedure in the realm of supersymmetric models should drive us to explicit 
higher order calculations to compute corrections in the MSSM, which may allow us to use 
our semi-nonpeturbative results in connection with the constraints on SUSY as imposed by 
the phenomenology of the LHC/ATLAS and CMS collaborations. 

The level of convergence of our results is satisfactory and, since our category of SUSY 
spontaneously broken model is still protected by the SUSY nonrenormalization theorem (the 
spontaneous breaking is a soft mechanism), we do not run into troubles with our perturbative 
calculations, for our coupling constants do not risk to take us to the strong coupling regime. 
This very same point must be reconsidered if we are dealing with a supersymmetric gauge 
theory, as it is the case of the MSSM, which is the ultimate framework to connect SUSY 
with the accessible energies. Before going directly to the MSSM, we intend to extend the 
calculations we have performed here to include the Fayet-Iliopoulos D-term [2| models, in 
which SUSY is broken and gauge symmetry may also be. In this case, there appears a 
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number of non-trivial aspects in connection with superspace and supergraph calculations, 
such as the gauge choice (unitary or 't Hooft's gauge choices in superspace) and a rich 
structure of ^-dependent terms [36J in the sector of gauge superfield propagators [37(. The 
consideration of the Fayet-Iliopoulos models is clearly mandatory as a step prior to extend 
our analysis to the physics of the MSSM. We shall be next focusing on this specific step of 
our project. 

Physics in three-dimensional space-time has been acquiring special interest, in view of 
a very rich diversity of lower- dimensional Condensed Matter (CM) systems that can be 
approached by quantum field-theoretic methods. On the top of that, more recently, renor- 
malizable and unitary massive gravity models has driven the attention to planar gravity 
models. We know that SUSY may be connected to both types of systems, CM and gravity. 
Now, the realization of M = 1 SUSY breaking in three space-time dimensions is very special, 
since the structure underneath it is real and not complex (as it happens in M = 1, D = 4 
SUSY or M = 2, D = 3 SUSY). Also, renormalizability allows a higher-power scalar po- 
tential in D = 3, so that F-term SUSY breaking demands reassessment and a number of 
peculiarities show up. Also, a Fayet-Iliopoulos term is not just like in M = 1, D = 4 or 
Af = 2, D = 3 SUSY, for the gauge potential multiplet is spinorial in simple 3D-SUSY. So, 
in view of these nontrivial aspects and the potentialities of supersymmetric planar systems, 
we believe that there should be some interest in re-analyzing the methods and re-discussing 
the results of the present work. 
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Appendix: Superspace Integrals 

In this appendix, the superspace integrals appearing in section III are listed. 
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The integrals appearing in the one-loop diagrams of 0(S 2 ) are: 

J ± (9,9) = J ^ 12 $ (l)$o(2) [Dl{k)e\D\{k)5t 2 ] [D 2 2 (-k)9 2 D 2 (-k)St 2 ] 
= (16) W ; 

J 2 (9,9) = J d% 2 %(l) [Dl{k)e\Dl{k)5t 2 ] [9 2 2 D 2 2 (-k)9 2 2 D 2 2 (-k)5 4 12 ] 
= (16) 2 (F ) ; 

J 3 (9,9) = J d% 2 [e 2 D 2 1 (k)6 2 1 D 2 (k)St 2 ] [9 2 2 D 2 2 {-k)9 2 2 D 2 2 {-k)5\ 2 ] 
= (16) 2 ; 

M9, 9) = J rf 4 ^ 12 $o(l)$ 1 (2) [Dl{k)elDl{k)5t 2 ] [D 2 2 {-k)D 2 2 {-k)5t 2 ] 
= -(16) 2 k 2 (F )M ; 

J 5 (9,9) = J d% 2 %(l)$ 1 (2) [Dl{k)9lDl{k)8i 2 ] [D 2 2 {-k)9 2 2 9 2 2 D 2 2 {-k)5t 2 ] 
= (16) 2 (F ><^> ; 

J 6 (9,9) = J d% 2 %(l) [Dl{k)9lDl{k)5 A 12 ] [D 2 2 (-k)D 2 2 (-k)5t 2 ] 
= -(16)V(F > ; 

J 7 (9,9) = J d% 2 ^(l) [Dl{k)9lDl{k)5t 2 ] [9 2 2 Dl{-k)Dl{-k)5t 2 ] 
= -(16) 2 k 2 ( Vl ) ; 

J 8 (0,9) = J d% 2 ^{\) [Dl{k)9lDl{k)8{ 2 ] [9 2 2 Dl{-k)Dl{-k)9 2 2 9 2 2 5t 2 ] 
= (I6) 2 (^i) ; 

J 9 (9,9) = J d% 2 [9\Dl{k)9lDl{k)5t 2 ] [D 2 2 (-k)D 2 2 (-k)St 2 ] 
= -(16) 2 fc 2 ; 
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Jio(M) = / ^^(1)^(2) [Dl{k)elD\{k)5t 2 ] [D 2 2 {-k)D 2 2 {-k)9 2 2 8\ 2 ] 

-(16) 2 £;%i> 2 ; (108) 



J n (e,9) = J d 4 ^ 12 $ 1 (l)$ 1 (2) [Dl{k)Dl{k)5\ 2 ] [D 2 2 {-k)D 2 2 {-k)5 A 12 ] 

= ; (109) 



Jn(0,9) = J d 4 ^ 12 $i(l)<f 1 (2) [Dl{k)Dl{k)6lel5t 2 ] [D 2 2 (-k)D 2 2 (-k)9 2 2 6 2 2 5l 2 ] 

= (16)V> 2 ; (no) 



J l3 (9,9) = j d% 2 ^(l) [Dl(k)Dl(k)5t 2 ] [Dl{-k)Dl{-k)5t 2 ] 

= 0; (111) 



J 14 (9,9) = J d% 2 [D 2 (k)D 2 (k)5$ 2 ] [D 2 {-k)D 2 2 {-k)5{ 2 ] 

= ; (112) 

Ji 5 (e,e) = J d 4 £ 12 $ (l)<5 (2) [Dl{k)5t 2 ] [D 2 2 (-k)5t 2 ] 

= 16(F ) 2 ; (113) 

Jie(M) = J ^i 2 $o(l)$o(2) [D 2 (k)5 4 12 ] [D 2 2 (-k)D 2 2 (-k)mD 2 2 (-k)5t 2 ] 

= (16) W ; (114) 

J 17 (9,9) = J ^ 12 $ (l)$o(2) [^(A;)/) 2 ^)^ 2 ^)^] [D 2 2 {-k)D 2 2 {-k)mD 2 2 {-k)5t 2 ] 
= (16) W ; (115) 

J 18 (9,9) = J d% 2 %{\) [D\{k)5t 2 ] [9lD 2 2 {-k)5t 2 ] 

= 16(F ) ; (116) 
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J 19 (9,9) = j d% 2 $ (l) [Dl{k)Dl{k)9lelD{{k)8\ 2 ] [9 2 2 D 2 2 (-k)D 2 2 (-k)mD 2 2 (-k)8t 2 ] 

= (16) 3 (F ) ; (117) 



JnfrS) = y d% [0?5?(fc)tf 2 ] [e 2 2 D 2 2 (-k)5t 2 ] 

= 16 ; (118) 



= (16) 3 ; (119) 



J2 2 (M) = y ^i 2 $o(l)$i(2) [Dl{-k)Dl{-k)Dl{-k)5t 2 ] 

= -(l6) 2 k 2 (F )( Vl ) ; (120) 



J 23 (9,9) = y d^ 12 $ (l)<l 1 (2) [D?(fc)£^(fc)^^(fc)^ 2 ] [Dl{-k)Dl{-k)Dl{-k)9l5 A 12 ] 

= -(16fk 2 (F )(^) ; (121) 

J 24 (M) = j d% 2 %{\) [Dli^Dli^lD^St,] [D 2 2 {-k)D 2 2 {-k)Dl{-k)9 2 2 5t 2 ] 

= -(16fk 2 (F ) ; (122) 

J 25 (9,9) = J d% 2 ^(l) [Dl{k)5\ 2 ] [9 2 2 D 2 2 (-k)9 2 2 D 2 2 (-k)D 2 2 (-k)5t 2 ] 

= -(Wfk 2 ^) ; (123) 

J 26 (M) = J d% 2 ^(l) [D 2 (k)D 2 (k)mD 2 (k)8l 2 ] [9 2 2 D 2 2 {-k)^D 2 2 {-k)D 2 2 {-k)5t 2 ] 

= -(16) 3 ^i> ; (124) 

J 27 (9,9) = J d% 2 [9 2 D 2 (k)D 2 (k)9 2 9 2 D 2 (k)5f 2 ] [ D 2 2 (-k)D 2 2 (-k)D 2 2 (-k)9 2 2 St 2 ] 

= -(lQfk 2 ; (125) 
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J 28 (M) = J d% 2 $ 1 (l)$ 1 (2) [Dl(k)6t 2 ] [D 2 2 (-k)5t 2 ] 

= ; 

J«(e,S) = I ^ 12 $ 1 (1)<I> 1 (2) [D\{k)5 A 12 ] [D 2 2 (-k)6 2 2 6 2 2 5t 2 ] 

= I6(^i) 2 ; 

J 30 (9,9) = J d% 2 ^(l)^(2) [Dl{k)Dl{k)9l9lDl{k)5t 2 ] [D 2 2 (-k)5t 2 ] 

= -(16) 2 k 2 ( Vl ) 2 ; 

J 31 (9,9) = J d% 2 ^(l)^(2) [Dl{k)Dl{k)9l^Dl{k)5t 2 ] [D 2 2 {-k)9 2 2 9 2 2 ^] 

= (16) 2 (^) 2 ; 

J 32 (9,9) = J d% 2 ^(l) [Dl{k)5\ 2 ] [D 2 2 (-k)6 2 2 6 2 2 5t 2 ] 

= 16(y?i> ; 

J 33 (9,9) = J d% 2 [D\{k)5t 2 ] [D 2 2 {-k)5t 2 ] 

J u (9,9) = J d% 2 [Dl{k)5\ 2 ] [D 2 2 (-k)9 2 2 9y i2 ] 

= 16; 

J 35 (9,9) = J d% 2 [Dl{k)Dl{k)9l9lDl{k)5{ 2 ] [D 2 2 (-k)5 4 12 ] 

= -(16) 2 k 2 ; 

J 36 (9,9) = J d% 2 [Dl{k)Dl{k)9l9lDl{k)5{ 2 ] [D 2 2 (~k)9 2 2 9 2 2 ^] 

= (16) 2 ; 
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J 37 (9,9) = J d% 2 $ 1 (l)$ 1 (2) [Dl{k)e\Dl{k)Dl{k)5t 2 ] [D 2 2 (-k)9 2 2 D 2 2 (-k)D 2 2 (-k)5 4 12 ] 

= (16) 3 A; 4 (^) 2 ; (135) 



J 38 (9,9) = j d% 2 ^(l) lDl{k)9lDl{k)Dl{k)5t 2 ] [D 2 2 (-k)6 2 2 D 2 2 (-k)D 2 2 (-k)5t 2 ] 

= W 3 k 4 M ; (136) 



J 39 (9,9) = j d% 2 [Dl{k)9lDl{k)Dl{k)5t 2 ] [D 2 2 (-k)9 2 2 D 2 2 (-k)D 2 2 (-k)St 2 ] 

= (16) 3 A; 4 ; (137) 



JaoWJ) = J d% 2 ^(l) [Dl{k)5\ 2 ] [D 2 2 (-k)5t 2 ] 

= . (138) 

The integrals appearing in the two-loop diagrams of 0(5 2 ) are: 



1,(9,9) = jd% 2 [Dl{p)Dl{p)9l5i 2 ] [D 2 2 (k)9 2 2 D 2 2 (k)5t 2 ] [^(^(-^fo)^] 

= 4(16)V ; (139) 

I 2 (9,9) = Jd% 2 [Dl{p)Dl{p)5 A 12 ] [D 2 2 {k)D 2 2 {k)5t 2 ] [^(g)^-^ 2 /) 2 ^)^] 

= ; (140) 

1 3 (9,9) = Jd% 2 [Dl(p)Dl(p)5$ 2 ] [D 2 2 (k)9 2 2 9 2 2 D 2 2 (k)5t 2 ] [^^ 2 2 (-#^(<?)^ 2 ] 

= 4(16)V ; (141) 

1 4 (9,9) = Jd% 2 [Dl{p)9l9lDl{p)8t 2 ] [D 2 2 (k)9 2 2 9 2 2 D 2 2 (k)5t 2 ] [^(^(-g)^ 2 ^ 4 ,] 

= -4(16) 3 ; (142) 

1,(9,9) = Jd% 2 [Dl{p)9lDl{p)Dl{p)8t 2 ] [D 2 2 (k)5t 2 ] [D 2 (^)D 2 2 (- g )D 2 (g)D 2 (g)^5 4 2 ] 

= (16)W ; (143) 
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X 6 (M) = Jd% 2 [Dl(p)6lDl(p)Dl(p)5t 2 ] [Dl(k)Dl(k)mD 2 2 (k)5t 2 ] 
x [m)Dl{-q)Dl{q)Dl{q)el5t 2 ] 
= (16) 5 pV ; (144) 



l 7 (6,e) = jd% 2 [Dl(p)5$ 2 ] [D 2 2 (k)5t 2 ] [Dl{q)Dl(-q)6l6l5t 2 ] 

= (16) 2 ; (145) 



U6,6) = jd% 2 [Dl(p)5 4 l2 \ [Dl{k)Dl{k)e 2 2 elDl{k)5t 2 ] [5 2 (g)£ 2 V^i W 2 ] 

= (16) 3 ; (146) 



1 9 (6,6) = jd% 2 [Dl{p)Dl{p)6l9lDl(p)8i 2 ] [Dl(k)Dl(k)mD 2 2 (k)5t 2 ] 
x [Dl( q )Dl(- q m5l 2 ] 
= (16) 4 ; (147) 

2x0(0,0) = Jd% 2 [Dl{p)8{ 2 ] [Dl{k)5{ 2 ] [Dl{q)Dl{-q)5t 2 ] 

= ; (148) 

X 11 (9,9) = Jd% 2 [Dl{p)5{ 2 ] [Dl{k)Dl{k)9l9lDl{k)8{ 2 ] [D 2 (q)D 2 2 (-q)St 2 } 

= -(16) 3 A; 2 ; (149) 

2 12 (M) = jd% 2 [D{{p)Dl{p)e{elDl{p)8i 2 ] [Dl{k)Dl{k)eldlD 2 2 {k)5t 2 ] 
x [Dl(q)Dl(-q)5t 2 ] 

= -(16)V . (150) 
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